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INCLUSION OF AERODYNAMIC 
DRAG IN VINTI'S SATELLITE THEORY 
Stan Watson 
N. L. Bonavito 

ABSTRACT 

In order to achieve separability in the Vinti theory of earth satellite motion 
when a non-conservative force such as air drag is considered, it is presently 
necessary to introduce a set of variational equations for the orbital elements, 
expressed as functions of the orbital elements and the tangential, radial, and 
normal components of the non- conservative forces acting on the system. The 
presence of a velocity dependent potential no longer allows the Hamiltonian to 
be the total energy, and furthermore, this term is cubic in the velocity. In this 
approach, the Hamiltonian is preserved in form, and remains as the total energy 
but the initial conditions and hence the Jacobi constants of the motion are changed 
with time through the variational equations . Hence separability is achieved at 
each point in time corresponding to continuously changing generalized coordinates 
in phase space. 

Some initial results of application of this program to the heavy air drag satellite 
San Marco 2, show excellent improvement in the accuracy of the computed orbits. 
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INCLUSION OF AERODYNAMIC 


DRAG IN VINTI’ S SATELLITE THEORY 


INTRODUCTION 

In this paper we shall treat the variation, of Izsak elements caused by dis- 
turbing forces acting on an orbit. These elements, an intrinsic part of the Vinti 
theory, are obtained as functions of time and then by a process of re-initialization, 
the Vinti equations of motion are solved for the position and velocity of the satel- 
lite, We should note that the air resistance of a body has in general six com- 
ponents, three being forces and three being moments of forces, which tend to 
make the motion of an asymmetric body very complex. However, for the case of 
a non-rotating sphere, the resistance can be reduced to a single component di- 
rected oppositely to the velocity of the sphere. In this treatment we shall intro- 
duce semi-empirically determined components of drag force which are functions 
of the eccentric anomaly. With use of these, the differential equations for the 
variation with time of the Izsak elements can be solved along with the Vinti equa- 
tions of motion. 

I. DRAG RATE EQUATIONS 

Here we present equations for the time rate of change of the Izsak orbital 
elements of tbe Vinti satellite theory (Reference 2). Since they can be shown to 
result from a consideration of conservation of energy and angular momentum, 
these equations hold for all perturbing forces. In our treatment, we shall assume 
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that these equations represent instantaneous departures due to .aerodynamic drag 
or perturbing forces on a Vinti orbit. These departures in the Vinti-Izsak orbital 
elements are then integrated as functions of eccentric anomaly for the instant of 
time under consideration. 

The differential equations representing variations of orbital elements are : 


da 

_2 dU 



dt 

na dM 



de 

1 -e 2 dU Vi- e 

l H 


dt 

na 2 e ^ na 2 e 

dco 


dco 

- cos i i 


YT ^ 

dt 

na 2 il - e 2 sin l 

vdl J + 

na 2 e 

di 

COS 1 / 

'*U\ 


dt ” 

na 2 Vi- e 2 sin i 

ictcoj 

/ 

dQ 

1 ( 

r du \ 


dt ~~ 

na 2 /I - e 2 sin l ' 

■ di 1 


dM 

1 ~e 2 (dtA 

2 dU 


dt 

n na 2 e 

na da 



( 1 ) 


where n is related to a by n 2 a 3 = GM = /a. G is the gravitational constant, and 
M is the mass of the Earth, u, is any disturbing function due to the oblateness 
of the Earth. In order to get these equations in terms of a perturbing force F, 
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we assume that this perturbing force may be described by the gradient of a 
perturbing potential. We now resolve the perturbing force F into the following 
components: 

R, in the direction of the position vector from the force center to the satellite. 

T, is perpendicular to R, lies in the orbital plane, and is positive in the di- 
rection of motion. 


W, is mutually perpendicular to R and T and completes a right handed bet of 
component directions. 

The partial derivatives appearing in Equations (1) are given by: 


£U 

d a 



dU 

d e 


- R a cos v + T 


2 + e cos v \ 
1- e 2 / 


r sin v 


dU 

da 


r W sm i p 


dU 

do. 


T r cos l - W( r cos sm i) 


(5U 

dco 


r T 


d U 



sin v + T 



( 2 ) 
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where \ p is the argument of latitude. Substituting (2) into (1) 


da _ /a 3 \ 1/2 1 r , 

3t " 2 \JJ >- Re sin v + T( 1 + e cos v)J 


^ /l - e 2 [r sin v + T(cos v + cos E)] 


r cosi/i 


di _ 

dt ” W \na 2 7T^ 


d co 
dt 


At- ~ R 


y 1 “ e 2 cos v 


n a e 


Vl - e 2 (l + — j sinv 


+ T 


n a e 


/ r cos i sin p \ 
Ina 2 /I - e 2 sini/ 


dn 

dF = " 


r sin p 


na 2 yi- e 2 sin i 


dM 

dt 


4l-e 2 ) 


2r 

cos v r | - T 

na / 


(l - e 2 ) ^ 1 + — ^ sin v 


nae 


(3) 


The Vinti parameter related to the orbit inclination is 


S - sm 2 i 


(4) 


From this, we have 


dS _ di 

dt “ 2 sm i cos i ' 
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or, 


dS 

dt 


__ di 


( 5 ) 


The variation of the inclination can then be written as 


dS 

dt 


2 /S( 1 ~ S) r\V cos \Jj 
//7a /l - e 2 


(6) 


Let Q - fi 3 , & = fi 2 , and using 


P - 



n 2 a 3 = fj. , 


sm i 


v§, 


we have, 


d h 

dt 


f/S [aR(l -e 2 ) cos v - r(2 + e cos v) T sin vj + rWe sin i/'/i - s"| 


(/i a) 172 eVS /l ~ e 2 


and 


d A 

dt 


r W sin »/> 

(/ia) 1/2 /S /l * e 2 


( 7 ) 
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The mean anomaly is related to the time of perigee passage /? r by the expression. 


M - n (t + /3 X ) 


( 8 ) 


From this. 


or, 


dM _ dn , , d ^i 

dF “ dt + ?i) + n ~dF + n 


d/3 1 

“dF 


1 f dM , s dn 1 

" n LdF ” ( fc H A) dt “ n J 


(9) 


and 


— - _ . 5 . .. 1/2 5/2 ^2 

dt 2 ^ a dt 


(10) 


together with dM/dt from (3), we find 


d/Sj 

dt 


da 


1 a ~ 1 dt + M(~2a 1 )~ 3/2 t fj. 


,1/2 o-5/2 „ _ 

3 dt 


da 2r R 


O a ) 


1/2 


(l-e 2 ) t*\ in 




R cos v - T 1 + 


(i 


a l - e 


’ (11) 


Summarizing our equations for the time rate of change of orbital elements, we 
have: 


da 

dt 


= 2 -7 


/ a 3 \ 1/2 l 

(jrj -^y==-[Re sinv +T(1 +e cos v)] 
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de 

dt 


/ a \l/2 

= T 1 - e 2 j^R sin v +T(cos v + cos E)J 


dS 

dt 


2y |f S( 1 ~ S) rW cos <p 
(/u,a) 1/2 Vl - e 2 


dt 


rW sin<p 


(/ia) 1/2 /l~e 2 


d^ 2 

dt 


V§ [aR( 1 “ e 2 ) cos v - r(2 + e cos v) T sin v] + reW sin ip Vl S 

1 A /r* ’ /" IT 


(^) 1/2 yge/l - e 2 




d/3, 

dt 


da 


da 2r R 


_ — B a -l + u.(~2a 'l"3/2.JdL t 1/2 a - S/2 , _ 

2 a dt + Aa XI ]2 a dt / 1/2 


(A^) J 


(l-e 2 ) 


(?) 


1/2 


R cos v - T [ 1 t 


a (l - e J 


( 12 ) 


n. SOLUTION OF DRAG-RATE EQUATIONS 

Let us now substitute for R, T, and W in terms of the eccentric anomaly as 
given m Reference 1, page 165* 


1 A dE 

R = ~ 2 C d¥ ae ^ V sinE dt 


1 A 


T = ~ 2 C D -(l-e 2 )^a P V 


1-d 


(1 - e cos E) 2 
(1-e 2 ) 


dE 

dt 


( 1 ) 
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w 


1 A 


1/2 a 3/2 


= ~ 2 C D¥ a ^s Y M a 


x ( 1 — e cos E) 2 l^s cos i jj 


dE 

dt 


where V is the velocity of the satellite defined by the equation 


1/2 

V = |-^-j (1 +e cosE)~ 1/2 (1 ~e cosE)“ 1/2 j(l +e cosE) 


- d( 1 + e cos E)J (2) 

where 

6 = co s jjT 1 ' 2 a 3 -' 2 (1 - e 2 ) 1/2 (1 - S) 1/2 (3) 

where 

<x> s = angular velocity of rotation of earth 
p — atmospheric density 
A = projected area of satellite 
m = mass of satellite 
C D = drag parameter {usually around 2„2) 

cos E - e 

cosv = 1 - e cos E 
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sm v 


(1 - e 2 ) 1/2 sinE 

r-e cos E 


r - a(l — e cosE) . (4) 

Now substituting Equations (2) and (3) into Equations (1), and inserting these 
results with Equations (4) into Equations (12) of Section I, we get 

-jj: = “ “ pa 2 -|( 1 - e cos E) _ ( 1 + e cos E) 1//2 [(1 + e cos E) 

, .-i dE 

“ d(l - e cos E)] 2 j dt 


= ” ^ ”rn P a (l “ e 2 ) ^ (1 + e cosE)“ 1/2 (1 - e cos E) 3/2 [(1 + e cos E) 


-- d( 1 - e cos E)] J + |e sm 2 E + ( 1 - e 2 ) 1 [( 1 - e 2 ) 


“ d(l - e cosE) 2 ] [(cos E~ e) + cos E(1 ~ e cosE)J ^ 


dE 

dt 
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dS 

dt 


= - Cq-p^ yS(T“S) a 5i/2 (1 -e 2 )' 172 


dE 


jcos 2 1 - e cos E) 5/2 ( 1 + e cos E)~ 1/>2 j^( 1 + e cos E) - d(l - e cos E)]| 


d/3 2 

dt 


1 A ( 

- ~ C D — pae~ l ( 1 “ e 2 ) -1/2 -je sin \ p cos </*(! ~ S) 


1/2 ^ /x~ 3/2 a : 3/2 

s ^ 


x (1 + e cos E)“ 1/2 (1 “ e cos E) 5/2 £( 1 + e cos E) 


- d(l - e cos E)] + (1 - e cos E) 3/2 (1 + e cosE) 1/2 [(1 + e cosE) 


- d(l - e cos E)j [2 sin E(1 - d) + de 2 sm E(1 + e 2 cos 2 E) +de sinE cos E 


1 dE 

4e 2 smE(l + d cos 2 E) * 2e 4 smE(5 - 2e 2 + e 2 cos 2 E)J j 


d (3 1 a 

= -=- r — m ,,-1/2 „5/2 , , _ p 2 r l/2 

dt 2 C D m f 1 a V 1 e ) 


x 


jsm i//cos >//(!- e cos £) 5/2 (1 + e cos E)~ 1/2 [(1 f e cos E) 


- d(l - e cos E 
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-M- 2a 1 )-V2^ t ^ 


+ (1 -e 2 ) 1/2 -JJT- + Vl -S 


A r* 

+ C D ~ pae lsinE(l - e cos E) 1/2 (1 + e cosE) - 1/2 [(1 + e cos E) 

->1 cIe! 3 1. „ . da 

“ d(l - e cos E)]J dFp 2 a (% "A) dt ‘ 


We now integrate Equations (5) from time t to time and indicate the 
eccentric anomaly at time tj by E 1 and the eccentric anomaly at time t 2 by E, 
We have the following equations : 




p {(1 - e cos Ey 1/2 (1 + e cos E)~ 1/2 [(1 + e cos E) 


- d (1 - e cos E)] 2 } d E 


1 ^ A 


^ e — m (1 ~ e2 ^ p {[(1 + e cos E')" 1 ^ 2 (1 - e cos E) 3/2 [(1 -s- e cos E) 


- d (1 - e cos E)]] [e sm 2 E + (1 - e 2 ") 1 [(1 - e 2 ) 


d (1 - e cos E) 2 ] [(cos E - e) + cos E (1 - e cos E)]] } d E 
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AS = -C '^co S vT^S /T 1/2 a s/2 Cl - e 2 r 1/2 

u m s 


x J ( ocosi//{(l - cos E'jS/ 2 (1 e co sE) [{1 +e cosE)-d(l - e cos E)] } d E 

E. 


/•Ei 


A ^ 2 = -^ D - a 6 "^ 1 -® 2 )" 172 
^ 2 u m 


,o {e stn t/>cos 0(1 -S) 1 / 2 /x 3 ^ 2 a 3 / 2 


x (1 + e cos E)“ 1//2 (1 — e cos E) s/2 [(1 f e cos E) 


- d (X - e cos E)] j- (1 - e cos E")~ 3//2 (1 + e cos E)“ 1,/2 [(1 + e cos E) 


- d (1 - e cos E)] [2 sin E (1 - d)- 4 a 2 sin E (1 f- d cos 2 E) + 2 e 4 SJLn E 


+ de 2 sin E {1 + e 2 co s 2 E) +de sin E cos E (5 - 2 e 2 i- e 2 cos 2 E)] } d E 


A/3 =i rAe, fT 172 a 572 (1 - e 2 r 172 
72 2 u m *“ 


r e 2 


p sin i/< cos i/; (1 - e cos E) 5 ^ 2 (1 ^ ecosE) 1/2 [(1+ecosE) 


-d(l-ecosE')] dE 
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A jSj = - P (- 2 c^y 3 / 2 j| t /^/2 a -s/2 A 


4 - (1 - e 2 ) 1 / 2 [A/3 2 + (2 - v/ 1 - S) A j0 3 ] 


+ v, - a e . 
u m 


p sin E (1 - e cos E) 1/2 (1 + e cos E)” 1/2 [(1 t e cos E) 


- d (1 - e cos E)3 d EV i (t 0 - £ ) A 

J jL 


( 6 ) 


Since p is a function of the eccentric anomaly also, then equations (6) are of the 
form 


p (E) d G (E), 


which can be integrated by parts. Following this, we can then write the solutions 
of equations (6) as, 


I 

e 2 • 

P (E) d G (E) = p (E) G (E) 

- 


E, 'E 


G (E> d p (E) 


(6. a) 


In order to find the functions G(E), we expand equations (6), keeping 0 (e) only as 
a factor. From this we obtain the terms to be integrated, 


(1 - e cos Er 1/2 (1 + e cos K ) s/2 d F 


(?) 
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(8) 



(1 " e cosE) -1/2 


( 1 + e cos E ) 1/2 dE 


(1 -e cos E) 1/2 (1 +e cosE)" 1/2 dE 


( 9 ) 


r*. 


J E, 


( 1 - e cos E) 1/2 ( 1 +■ e cos E) 3/2 dE 


( 10 ) 


r 

Jr 


(1 -e cosE ) 3/2 (1 +e cosE) 1/2 dE 


( 11 ) 


(1 “e cos E) 5/2 (1 + e cosE)' 1/2 dE 


( 12 ) 


(1 - e cos E) 5/2 (1 +e cosE) 1/2 dE 


( 13 ) 


J E, 






( 14 ) 


f 2 

cos E( 1 - e cos E) 5/2 ( 1 + e cos E) 1/2 dE 


( 15 ) 
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(16) 



sin E(1 - e cos E) s,/2 


( 1 + e cos E)~ l ' 2 dE 


rE. 


J E . 


sm E( 1 - e cos E)~ 3/2 ( 1 + e cos E ) 1/2 dE 


1/2 


(17) 



sin E( 1 - e cos E) 1/2 ( 1 + e cos E) 1/2 dE 


sin E(1 - e cos E) 3/2 (1 + e cos E)~ 1/2 dE 


sm E cos 2 E( 1 - e cos E) 5/2 ( 1 + e cos E) 1/2 dE 


sm E cos 2 E( 1 ~ e cos E)~ 3/2 ( 1 + e cos E) 1/2 dE 


(18) 


(19) 


( 20 ) 


( 21 ) 


r E 2 

I sm E cos E( 1 ~ e cos E)~ s/2 ( 1 + e cos E) - 1/2 dE 


(22) 


r e„ 


sin E cos E( 1 ~ e cos E) 3/2 (1 + e cos E) 1/2 dE 


(23) 
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i ( 24 ) 



sin Ecos 3 E(1 - e cos E)" s/2 (1 + e cos E) -1/2 dE 


sm E cos 3 E( 1 — e cos E) _ 3//2 ( 1 + e cos E) dE 


cos E(1 - e cos E) 7/2 ( 1 + e cos E)“ 1/2 dE 


cos 2 E( 1 - e 2 cos 2 E) 1/2 dE 


cos E( 1 - e 2 cos 2 E) 1/2 dE 


( 25 ) 


( 26 ) 


( 27 ) 


( 28 ) 


f E2 

cos E sm E( 1 - e 2 cos 2 E) 1/2 dE 


( 29 ) 



sin E(1 ~ e 2 cos 2 E) 1/2 dE 


sm 2 E(1 - e 2 cos 2 E) 1/2 dE 


( 30 ) 


( 31 ) 
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(32) 


r E„ 


cos 2 E( 1 - e cos E ) 3 ' 2 ( 1 + e cos E)“ 1/2 dE 



cos E( 1 “ e cos E) 3/2 


(1 + e cos E) 1/2 dE 


(33) 


r E, 


cos E sm E( 1 - e cos E) 3/2 (1 + e cos E )~ l/ ' 2 dE 


(34) 


sm E(1 - e cos E) 3/2 (1 + cos E) 1/2 dE 


- 1/2 


(35) 


r E, 


'E, 


sin 2 E( 1 - e cos E) 3//2 (I + e cos E) 


V- 1/2 


dE 


(36) 


Some of the integrals (Equations 7-36) may be expanded into two or more 
integrals. Equation (7) may be used as an example. 
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(1 -e cos E)“ 1/2 (1 + e cos E) 5/2 dE 


2 

(1 “ e cos E) -1/2 (1 + e cos E) 4/2 (1 + e cosE) 1/2 dE 


( 1 - e cos E) 1/2 (1 4- e cos E) 1/2 (1 + 2e cos E + e 2 cos 2 E) dE 


( 1 - e cos E) 1/2 (1 + e cos JE) ly/2 dE 


r E 2 

+ 2 e cos E( 1 - e cos E)~ 1/2 ( 1 + e cos E) 1/2 dE 


e 2 cos 2 E( 1 -e cos E)" 1/2 (1 +e cos E) 1/2 dE.(37) 






18 



Using similar methods on Equations 8-36, we have the result 


(1 ~ e cos E) - 172 (1 + e cos E) 572 dE 


f E * 

j ( 1 - e cos E) 172 ( 1 + e cos E) 1/2 dE 


f E 2 

l-2j e cos E( 1 - e cos E)~ 1/2 ( 1 + e cos E) 1/2 dE 

- 


+ 



e 2 cos 2 E(1 -e cosE)~ 172 


(1 +e cos E) 1/2 dE (37) 


(1 - e cos E) 1/2 ( 1 + e cos E) 372 dE 



(1 - e 2 cos 2 E) 172 dE + f e cos E( 1 ~ e 2 cos 2 E) dE 


(38) 
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f E 2 

(1 - e cosE) 372 (1 +e cosE) 172 dE 


(1 -e 2 cos 2 E) 1/2 dE 


f 


cosE(l - e 2 cos 2 E) 1/2 dE (39) 



(1 - e cos E) 5/2 


(1 


+ e 


cos E)~ 1/2 dE 



(1 - e cos E) 172 (1 + e cosE) -172 dE 


- 2 



e cosE(l - e cosE) 1/2 


(1 + e cos E)' 172 dE 


f E 2 

e 2 cos 2 E(1 - e cos E) 1/2 (1 + e cos E)~ 1/2 dE (40) 

Jjf. - 
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- e cos E) 5/2 (1 + e cos E) 1/2 dE 


- r 


£ 


(1 - e 2 cos 2 E) 1/2 dE -2 e cosE(l - e 2 cos 2 E) 


, 1/2 


dE 


+ 


e 2 cos 2 E(l~ e 2 cos 2 E) V2 dE(41) 


(1 - e cosE) 7/2 (1 + e cosE) 1/2 dE 


■ r 


(1 -e cos E) 1/2 (1 + e cosE) 1/2 dE 


3 f e cos E(1 - e cos E) 1/2 (1 + e cosE) - 1/2 dE 

/e. 


+ 3 f e 2 cos 2 E(1 ~ e cos E) 1/2 ( 1 + e cos E) dE 
% 


f 


e 3 cos 3 E(1 - e cosE) 1/2 (1 + e cosE) 1/2 dE (42) 
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r E, 




cos E( 1 - e cos E) 5/2 ( 1 + e cos E) 1/2 dE 


1 _ 

e 


[r 


e cos E(1 - e 2 cos 2 E) 


, 1/2 


dE 


- 2 



cos 2 E(1 ~ e 2 cos 2 E) 1/2 dE 


+ 



cos 3 E(1 -e 2 


cos 2 E) 17 2 dE 


(43) 



sin E( 1 “ ecosE) 372 ( 1 + e cos E)" 1/2 dE 


sm E(1 “ e cos E) 1/2 ( I + e cos E) 


1/2 


dE 


f E 2 

e sin E cos E( 1 ~ e cos E) 1/2 ( 1 + e cos E) 1/2 dE (44) 
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f * 2 

sm Ecos 2 E(l-ecos E)“ 5/2 (1 + e cos E)~ 1/2 dE 



sm E(1 ~ e cos E) -2 ( 1 ~ e 2 cos 2 E)~ 1/2 dE 


/*E, 


; E, 


sin 3 E( 1 - e cos E)' 2 (1 -e 2 cos 2 E)~ 1/2 dE (45) 


sm Ecos 2 E(1 ~ e cos E)“ 3/2 (1 + e cos E) 1/2 dE 



sin E( 1 - e cos E) 3//2 


( 1 + e cos E) 1/2 dE 


s in 3 E( 1 - e cos E) 3/2 ( 1 + e cos E) 1/2 dE (46) 




sm Ecos 3 E(1 - e cos E)~ s/2 


( 1 + e cos E)~ 1/2 dE 


/*E, 


sin Ecos E(1 - e cos E)“ 2 (X - e 2 cos 2 E)~ 1/2 dE 


r 


sin 3 Ecos E(1 - e cos E) 2 ( 1 - e 2 cos 2 E)“ 1/2 dE (47) 
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sm Ecos 3 E(1 - e cos E )“ 3/2 


( 1 + e cos E) 1 ^ 2 dE 



s in E cos E( 1 “ e cos E) 3/2 ( 1 


+ e cos E ) 1/2 dE 


•'E, 


sm 3 E cos E( 1 - e cos E) 3/2 ( 1 + e cos E ) 1/2 dE 


(48) 


Below is the series solution for the incomplete elliptic integrals of the first 


and second kind. 

F(a, 4>) - 4> (l + J k 2 + 64 k 4 + -25Q k 6 + (64) (256) k8 ) 


1 

4 


k 2 sm cj> cos <fi { 1 + 


_ 9 _ 

16 


— ^4 

64 k 


49 

256 


25 

64 


1 / 25 , 

- 22 k4 sm 3 <j£) cos \3 + y^ k J 


1225 

12 



5 

96 


k 6 sin<£cos<£ 


( 49 

( 1+ 64 


2 - 


35 

"j 28 k8 sin? 0 cos (p (49) 


24 



E(a, 0) 


4> 


H 



175 

~ 16384 



+ ■ 32 ’ k 4 sin 3 0 cos 0 


+ 96 k 6 sm 5 0cos0 


(1 + vs k 2 + m k4 ) 

(-IH 


+ "-^Q 24 k 8 sin 7 0 cos 0 - (50) 

Here the values of k and 0 depend upon the term that is being integrated and 
these values can be found in Reference 3. In addition, throughout the report we 
make use of the following two substitutions in Equations (49) and (50): 

7 r 

0 - “ E + ~2 


and 


k 2 = e 2 
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Examples of how the integration is performed is given in the appendix. 
Here we give only the results of the integration. 


(1 + e cos E) 1/3 (1 - e cos E) 1/2 dE 


|ln [(I - e 2 cos 2 E) 1/2 + e smE] + F(a, <£)|. 


(51) 


f 

**' T? 


(1 - e cos E)~ 1/2 (1 + e cos E) s/2 dE 


In. [(1 - e 2 cos 2 E) 1/2 + e sinE] + F(a, <jb) 


t 2{ln [(1 ” e 2 cos 2 E) 1/2 + e sm e] 


+ E(a, <£)~F(a,4>)} + — In [(1 - e 2 cos 2 E) 1/2 + e sinE] 


e sinE(l - e 2 cos 2 E) 1/2 

2 +E(a, cp) - F(a, <£) 


7 €• ^ 

~2 In [(1 - e 2 cos 2 E) iy2 + e sinE] 


l oTp/ , N -i v., e smE(l - e 2 cos 2 E') 1/2 l 

+ 3 [E(a, 0) -F(a, <p)j F(a, <£) — 2 f 

J 


(52) 
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f 2 

l (1 “ e cos E) 1/2 (1 + e cos E)' 1/2 dE 
Je 1 


= {- 2F (a, <£) + E(a, <£)} 


(53) 


r E 2 

l ( 1 -e cos E) 1/2 (1 + e cos E) 3 2 dE 

JEi 


- E(a, <£) + sinE(l - e 2 cos 2 E) 


, 1 / 2 


+ (l-e 2 )ln 


e sinE + (l~e 2 cos 2 E) 1 /2 
(1 - e 2 ) 1/2 


(54) 



(1 - e cosE) 3/2 (1 


+ e cos E) 1-/2 dE 


-E(a, ft 
_ 


- \ |e sinE(l “ e 2 cos 2 E) 1/2 


+ 


(l-e 2 )ln jesinE + (l~e 2 



(55) 
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( 1 - e cos E) s/2 ( 1 + e cos E) 1/2 dE 


F(«, 4 >) + 


e sinE(l -e 2 cos 2 E) 

2 


7 + e 2 
2 


In [( 1 - e 2 cos 2 E) 1/2 + e s in e] 


+ 


3E(a, 4>)~ 3F(a,‘ 4>y. 



( 1 - e cos E)- 7/2 ( 1 + e cos E )~ 1/f2 dE 


F(a, 4>) - 2(3 + e 2 ) In [(1 - 


e 2 cos 2 


E) 


1/2 


+ e sin E] + 8E(a, (p) — 8F(a, <f>) 


+ 2e sinE(l 


~ e 2 


cos 


: E) 1/2 - 3 


— 2 3 

e cos 


Mi 


~ e 2 


cos 


sin E 




1/2 


_1 

3 



^,)-(l~e 2 )F(a, <P) 


(1 + e 2 ) 


sm 4> cos <fi 
/l - e 2 sin 2 <£ 


cos <$ Vl ~ e 2 sin 2 <p 


sm <fi 


(57) 
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cos E( 1 - e cos E) 5/2 ( 1 + e cos E) 1/2 dE 


5 +3e- 


e sinE(l " e 2 cos 2 E) 


2 ™ 0 2 pM /2 


+ (1 - e 2 ) In 


sinE+(l~e 2 cos 2 E) 1,/2 

(1 _ e 2 ) l /2 


'2 fl -e 2 cos 2 E') 1/2 + ecos 3 E 2j„ 

+ 3 JIHe +3-jE(i,0)-(l-e 2 )F(», 


2 (1 + e 2 ) sin0 cos 4> 

Vl - e 2 sin 2 4> 


2 ) s in cf> cos 4> cos Vl ~e 2 sin 2 <^> 


sin 0 


- x e sinE(l -e 2 cos 2 E) 3/2 


r E 2 

sin E( 1 - e cos E) 5/2 ( 1 + e cos E)“ 1/2 dE 

- 


3e( 1 **e 2 cos 2 E) 3/2 


1 2 

{2 +2e cosE +e cosE(l -e 2 cos 2 E)j (59) 

IT 
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smE(l-e cos E)“ 3/2 (1 + e cos E) 1/2 dE 


e(l - e 2 cos 2 E) 


1/2 


j^2(l + e cosE) 


+ (l~e 2 cos 2 E) 1/2 sin -1 [( 1 - e 2 cos 2 E) 1/2 ] } 


(60) 


sin E( 1 ~ e cos E) ly ^ 2 ( 1 + e cos E) 1//2 dE 


•'Ei 


= 9^- jsm 1 [( 1 - e 2 cos 2 E) 1/2 


- e cosE(l - e 2 cos 


E) >'»]} 


( 61 ) 



sin E( 1 ~ e cos E) 3/2 


(1 + e cos E) -1/2 dE 


r~ 


1 

e 



[( 1 - e 2 cos 2 E) 1/2 ] 


- 2 ( 1 ~ e 2 cos 2 E) 1/2 + ~2 e cos E (1 ~ e 2 cos 2 E) 


. 1/2 


(62) 
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r E a 

J sin E cos 2 E( 1 - e cos E)~ 5/2 (1+ e cos E)~ 1/2 dE 

E i 


3e(l -e 2 cos 2 E) 3/2 + 2ecos ® + e cos E(l- e 2 cos 2 E) 


e 2 [“ 3(1 - e 2 cos 2 E) 3/2 sin" 1 [(1 - e 2 cos 2 E) 1/2 ] + 


ecos E(l- e 2 ) (1 - e 


* C 1 + ecosE) [2(1- e 2 ) ~ 6(1 - e 2 cos 2 E)]]J 


(63) 


sm E cos 2 E( 1 - e cos E)" 3/2 ( 1 + e cos E) 1/2 dE 


e( 1 - e 2 cos 2 E) 1/2 | 2(1 + ecos E) 


+ (1 e 2 cos 2 E) 1/2 sin* 1 [(1 - e 2 cos 2 E) va ‘ 


2 gS" j(S -2e 2 ) sm -1 [(1 - e 2 cos 2 E) 1/2 ] 


4(1 -e 2 ) (1 + e cos E) ( 1 - e 2 cos 2 E) - 1/2 


+ 4(1 -e 2 cos 2 E) 1/2 + e cos E( 1 - e 2 cos 2 E) 1/2 


(64) 


2 cos 2 E) 
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sin E cos £( 1 - e cos E) ’ 5/2 ( 1 + e cos E)~ 1/2 dE 


■ 1/2 


1 

3e 2 (1 - e 2 cos 2 E) 3/2 


[3(1 - e 2 cos 2 E) - 2 - 2e 3 cos 3 e] 


E 


2 


(65) 


E 


1 



sin E cos E( 1 - e cos E)~ 3/2 


( 1 + e cos E) 1/2 dE 


e 2 (1 - e 2 cos 2 E) 


1/2 


■^2(1 + e cos E ) + ( 1 -• e 2 cos 2 E) 


+ 2(1 - e 2 cos 2 E) 1/2 sin -1 [(1 - e 2 cos 2 E) 1 ^ 2 ] | 


E 


2 


( 66 ) 



s in E cos 3 E( 1 - e cos E)“ 5/2 


(H 


e cos E) -1//2 dE 


1 

3e 4 (1 - e 2 cos 2 E) 3/2 


^8 + 8e cosE(l -e 2 cos 2 E) 


- 3e 2 cos 2 E( 5 ~ e 2 cos 2 E) - 2e 5 sin 2 E cos 3 E 


+ 6(1 " e 2 cos 2 E) 3/2 sin" 1 [(1 


e 2 cos 2 



(67) 
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sinEcos 3 E(l~ ecos E)~ 3/2 (1 + e cos E) 1/2 dE 


(1 


e 2 cos 2 


E) 


1/2 


|^2(1 + e cos E) 


+ ( 1 - e 2 cos 2 E) 


+ 2(1 -e 2 cos 2 E) 


1/2 sm 1 


[ci- 


e 2 cos 2 


E) 


I/2 1] 


i r 

^[ (1 - 


e 2 cos 2 E) 1/2 (3 - e 2 + e cosE) 


+ (3 - 2e 2 ) sin -1 [(1 - e 2 cos 2 E) 1/2 



(1 -e 2 cos 2 E) 3/2 


+ 2(l - e 2 )(l-e 2 cos 2 E) _1,/2 ( 1 + e cos E 
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cos E( 1 - e cos E) 7/2 


(1 + e cos E)~ l/2 


1_ 

e 


e sin E( 1 ~ e 2 


cos 2 E) 


3/2 


+ e sinE(l ~e 2 cos 2 E) 


1/2 


-ft (l-e 2 )-4 


3^ 

8 


(1 -e 2 ) 2 -4(1 -e 2 ) +8 


I ln[(l- 


0 2 
e cos 


E) 


1/2 


+ e s m 


E ] 


_4 (1 - e 2 cos 2 E) *^ 2 + e 2 cos 3 E 
3 sm E 


- 8E(a, <f>) + 8F(a, <j>) 



E(a, <P)~ ( 1 “ e 2 ) F(a, <f>) 


- 


_ simpcosd) cos<£(l- 

(1 + S } M -o 2 sm 2 ^ 2 + 


2 j \ 1/2 


(69) 



cos 2 E(1 _ e 2 


2 

cos -4 


E) 1/2 dE 


1 e 2 cos 3 E(1 ~ e 2 cos 2 E)^ 2 
e 2 i 3 sinE 


(1 - e 2 ) 
+ 3 



E 


2 


E 


1 


( 70 ) 
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cos E(1 - e 2 cos 2 E) 172 dE 




e sinE(l ~e 2 cos 2 E) 


1/2 


+ 


(1 - e 2 ) In [e sinE + (1 " e 2 



(71) 


f 

•'f*. 


2 1 
cos E sm E(1 - e 2 cos 2 E) 1/2 dE - 


3e 2 


— r'rtc 2 3/2 


( 1 - e 2 cos : 


(72) 


sm E(1 - e 2 cos 2 E) 172 dE 


- 2 e^^ sin " 1 [Vl - e 2 cos 2 eJ - e cos E Vl - e 2 cos 2 


(73) 


f E 2 

J sin 2 E (1 - e 2 cos 2 E) 172 d E 


i. ( 1 ~ e 2 cos 2 E) 1/2 sin 3 E 
3 cos E 


f 


E(a, 4 >~) + 


cos 4>1l -e 2 sin 2 


0 


sin 0 


~ \ [( 1 -e 2 )F(a, 0) “ E(a, 0)] 


( 74 ) 
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cos 2 E(1 - e cos E) 3/2 (1 + e cos E)~ 1/2 dE 


2 r 1 f r 

2 E(a, <£) - F(a, <fi) + 2 |e sinE Tl - 

L_ 


e 2 cos 2 E 


+ (1 - e 2 ) ln[esmE+ Vl - e 2 cos 2 e]T - In [e sin E + Vl - e 2 cos 2 e] — ^-©1 

J E i 


cos E( 1 - e cos E) 3/2 (1 + e cos E)~ 1/2 dE 


2e *^- n [ e Sin E + il ~ e 2 cos 2 e] ~E(a, <£) + F(a, cfi') 

- E 

lj~ jje s In E /l - e 2 cos 2 E + ( 1 - e 2 ) In [e s m E + i 1 - e 2 cos 2 Ejjj- (76) 


cos Esin E(1 - e cos E) 3/2 (1 + e cos E)~ 1/2 dE = ^2 Vl - e 2 cos 2 E 


e 2 cos 2 E J r (77) 

J E, 
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f 


s in E( 1 - e cos E) 3 ' 2 ( 1 + e cos E) - 1/2 dE - s ^- n ~ 1 “ e 2 cos 2 e] 


- 4 Yl~ e 2 cos 2 E*ecosEVl ~ e 2 cos 2 E^ 


( 78 ) 


rE, 


sin 2 E( 1 - e cos E) 3 ^ 2 ( 1 + e cos E) - i//2 dE 



e 2 ) F(a, 4>~) - E(a, <£)] 


+ (1 - e 2 ) In [e smE + Vl - e 2 cos 2 e] 


- e sin 


E - /! “ e 2 cos 2 E-©} 


(79) 


here 


© 


- E(a, <£) + F(a, <£) 


+ 


|E{a, <£) + e 2 F(a, <p) 


sin <£71 — e 2 sm 2 <£ 
cos <£ 
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Equations (51 through (79) give us the functions G(E) of equation (6 .a). For 


P( E), we adopt the expression given by King-Hele (Reference 4), 


p (E) - p (1 + b x 2 (1 - cos E) 2 } e 


mxl-cosE) 


(80) 


where p p and H p are the density and density scale height at perigee respectively, 
x = ae, and b = p./2H 2 [l +(l/2) pj. Here p is a constant which represents the 
rate of increase of H the density scale height, and is most likely to have a value 
near 0.1. The scale height is defined (Reference 5) as 

h 

H = E (81) 

o 

where R* is the universal gas constant, h p *is the pressure scale height given 
in the Tables, M 0 is the mean molecular weight of air at sea level, and the gravi- 
tational constant as a function of altitude is 


g = 


So 


r 2 


( r e tZ ) 2 


(82) 


Here is the effective earth's radius at a specific latitude, and Z is the geometric 
altitude. 

Now differentiating equation (80) with respect to the eccentric anomaly, we 
obtain the following result; 
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~(x/H p )( 1-cos E) 


sin E 


x 


c) P (E) = p p e 



+ 2 b x 2 - 


H 


~ 2 b x 2 cos E sin F [1 - 


H 


H 


co s 2 E sin E 


( 83 ) 


Equations (51) through (79) and (80) determine the first term of equation (6. a). 

We now use equation (83) together with (51) through (79) to obtain each integral 
of equation (6. a). For a fairly circular orbit, the integral term of (6. a) has a 
small effect since the atmospheric density does not change much. 

Multiplying each term of equation (83) by each of the results of equations 
(51) through (79), expanding each of the radical terms in a Taylor expansion such 
as, 


(lix)» = linx + n C n ~ x 2 

2 


(84) 


and 


(1 ± x)" n = 1 + 


n x + 


n (n -5- 




Expanding the term, 
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In (e sxn "5 -f- Vl - e 2 cos 2 E) 


= In [\Zl~T 2 ] + e 5111 S 

(1 - e 2 ) 1/2 


e 3 sin 3 E 


(1 _ e 2 ) 3 / 2 3 ! 


( 85 ) 


and retaining those terms in k 4 (e 4 ) m equations (49) and (50), we obtain the 
following set of integrals and solutions, 


r 


E sm Z d F. = sinE-Ecos E 


( 86 ) 


F 2 sm Ed E? 2 E sm F + 2 cosE 


- E 2 co s E 


(87) 


E 3 sin E d F = (3 E 2 - 6) sin E 


- (E 3 - 6 E) cos E 


( 88 ) 
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(89) 



E 4 sin E d E = {- E 4 co S E + 4 [(3 E 2 - 6) cos E 


+ E 3 sm E - 6 E sin 




E s sin E d E = 


{- E 5 cos E + 5 E 4 sm E 


(90) 


- 20 X (8S)}|* 2 


where I (88) denotes the result of integral (88) above 


* e 2 

sin F cos F d E =— sin 2 E 

J • 2 

E i 



(91) 



E sin E cos E d E = — F sin 2 


E 1 

£ - — + — sin E cos E 
4 4 


2 


E 2 sinEcos E = -E sin 2 E 
4 


--E 2 cos2 E 
4 



+ - cos 2 R 

8 


(92) 


(93) 


41 



( 94 ) 


E 3 sin E cos F d E = -1 e 2 sm2E-lsin2E 

8 16 


- If 3 cos 2 E + -E cos 2 E 
4 8 




E 4 sin E cos E d E = — [(4 E 5 - 6 E) sin 2 E 

8 


(95) 


- (2 E 4 - 6 E 2 + 3) cos 2 E] 



E 5 sin E cos E d E =— [(5 E 4 

8 


15 E 2 + 15 E) sin 2 E 


(96) 


- (2 E 5 - 10 E 3 + 15 E) cos 2E]| l 2 


sin 2 E cos EdE = _ sin 3 E 
3 


(97) 



sm E cos 2 E d E = - — cos 3 E 

3 


( 98 ) 
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sin E 



E 2 sin 2 EcosEdE= — EcosE+i- 

2 2 



(— - . 1 -1 s i n 3 E + cos 3 E 

\^12 54 j 18 


( 99 ) 



sin 2 E cos 2 Ed E = -A.E 

8 


— sin 4 E 
32 



( 100 ) 


f E 2 

E sin E cos 2 Ed F = I (86) 

J 

E i 


E cos 3 E 

12 


sin 3 E 
36 


3 E cos E 3 sin E 
4 * 

4 4 



( 101 ) 


I 


E 


2 


E 


1 


E 2 sm E cos 2 E d E = I (87) 




cos 3 E 


2e 2 f cos E f— E sin E (102) 

,42/2 


1 

18 


E sin 3 F 
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2 


F 4 sin 2 E cos E d E = I sin E '1 (95) 


_Il (103) + A I (101) + -£1(99) 
+ — sm E - ■§ I (97) -tEcosE 

5 5 5 


6 E 2 sin E + 


2e 4 sin E +1e 3 COS E 
5 S 


(105) 


IE, 


| 2 E s sin E cos 2 E d E = -^ cos 2 E I (90) 


+ 


10 

3 


I (105) - 40 I (99) + 80 I (97) 


+ ™ I (103) - 40 I (101) 
3 



(106) 


p B 2 

| E 4 sin E d E = 

% 


{ _ E 4 cos E + 4 [(3 E 2 - 6) cos E 


+ E 3 sin E 


- 6 E sm E]> . 


(107) 



sin 2 E d E =~ 


sin 2 E 
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(109) 


r* E q 




E 2 sin 2 E d E = — - 

-i) 

„ „ E cos 2 E 
sm 2 E - 

4, 

\ 4 

v 

4 


E 4 sm 2 E d E = ^ sin E I (107) + _ 2 E 3 


+ 61 (109) + 6 E + 3 sin 2 E 


- 2 1 (94) - 12 X (92) 


(ilO) 


sin 4 E d E = - E - I sm 2 E + A sm 4 E 
8 4 32 


( 111 ) 


r E 2 


sm E co s 3 EdE = -_ co s 4 E 

4 


l E , 


’E, 


( 112 ) 



E 2 sm E cos 3 EdE = 


3 

4 


cos E I (102) + 


125 

36 


I (91) 


2 

27 


I (112) +-L E 2 


— E sm 2 E - 
24 


— cos 2 E 
48 


-lEUUD-i 



(113) 


46 





E 4 Sin 2 E cos 2 E d E = E 4 I (100) 



— Le 3 cos 4E + -L.E cos 4 E - 
32 256 


3 

1024 


+ — E 2 sin 4 E 
128 


sin E cos 2 E \/l - e 2 co s 2 E d E = 

4 e 2 


e cos E 

2 


\Jl - e 2 ,cos 2 E 


- e 3 cos 3 E \/l - e 2 cos 2 E •+■ — arc sin (\/l - e 2 cos 2 E) 



E sin E sin 4 E d E = — E sin 3 E + ~ cos E 

3 75 


—cos 3 E - — E sm s E 

675 5 


O 

sin 4 E cos E 

25 


E 


2 


E 


1 



sin 4 E cos E d E = - 


cos 5 E 

~To 


cos 3 E 

6 


E 


2 


E 


1 




Ecos EdE=cosE + E smE 


sin 4 E 
(114) 


(115) 


(116) 


(117) 


(H8) 
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sin 2 Ecos 3 EdE = E cos El (100) 


-ll (118) + JL I (117) + 1 1 (87) 
8 32 8 


1 

32 


I (116) 



(119) 


E 2 sin E cos 4 E d E = 


cos 4 E 
~5 ' 


I (87) 



2 I (119) - 


2 cos 5 E 
5 


E 2 

E 1 


( 120 ) 


f E 2 


E 2 sin E \/l - e 2 cos 2 E d E = I (87) - — I (102) 

2 


-j I O20) 


( 121 ) 


F S 

E cos 6 E d E = — cos s E sin E + — E cos 3 E sin E 
6 24 


t rf c 1 

+ 1 — + ■ — E sin E cos E + — cos 6 E 
32 16 36 


+ <nr cos4E "Tif 1(911 


( 122 ) 


48 




s in E cos 5 E d E - 


E 2 cos 6 E 
6 


+ — I (122) 
3 


E i 


( 123 ) 


r 


E 2 sin E cos E y/ 1 - e 2 cos 2 E d E = I (93) 


- — I (113) -^-1 (123) 
2 8 


(124) 


r 


E cos 7 EdE-- cos 6 E sin E + ^ cos 4 E sin E 
7 35 


a 8 

j . cos 2 E sin E + — — sin E 

15 15 


E, 


(125) 


~ E ~ 


E 2 sin E cos 6 E d E = 


E^_co£_E 2 : (125) 
7 7 


(126) 


f 


E 2 sin E cos 2 E \J 1 - e 2 cos 2 E d E = I (102) 


— I (120) -i- I (126) 

2 o 


(127) 


49 




E3 cos 5 EdE = i- cos 4 E sin E 


+ — E 3 sin Ecos 2 E + ^-E 3 sin E 
15 15 


il (120) -il (102) 

5 ^ 


I (87) 





E 4 sin E cos 4 E d E = - 


E 4 cos 5 E 
5 


+ il (128) 


e 2 

El 


f 

•'IT 


E 4 sin E \/l - e 2 cos 2 E d E = I (107) 


e 2 e 4 


I (104) I (129) 

2 8 


(128) 


(129) 


(130) 
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2 F . 4 cos 6 E .E 3 sinEcos 5 E 

1 E 4 sin E cos 5 E d E = — ^ + “ g 

4, 


, 5 £4 , JL E 3 sm E cos E +-^ E 3 sin E cos 
+ 96 24 36 


- — I (123) (93) - ^ I (11^ 

3 8 12 


e 2 



E 4 cos 4 E 

E 4 sin E cos 3 E d E = - 7 4 



+ ^.F 3 sin E cos E + ~ sin E cos 3 E 
8 


| E 2 

_9 1 (93) _ 3 j (113) [ 

8 f El 


E 4 sin E cos E \fl - e* cost's d E = I (95) 

-£l (132) I (13D 
2 8 


p E 2 sin 3 E d E = - E 2 cosE+Ie 2 cos 3 E + |cos E 


+ lEsmE-ig Esin 3E-1- cos 3E 
2 


l E - 


E 


(131) 


(132) 


(133) 


(134) 
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_ _-> Jcos 6 E sin E 6 4 . 

E 3 cos 7 E d E = E J i + — cos 4 E sin E 

l 7 35 


+ i SinE -3^ sin3E f-f I(126) 


|| I (120) -|| I (87) + 1| I (134) 


(135) 


E 4 sin E cos 6 Ed E s - — cos 7 E + (135) 


(136) 


r 


E 4 sin E cos 2 E \Zl_e 2 cos 2 E d E = I (104) 


- ^ I (129) - ^ I (136) 

A O 


(137) 


f 


cos 2 E s/l - e 2 cos 2 E d E = 


- (1 - e2) {.» 
V 3 e2 1 


Si-ri ^ \(l ~ e 2 sm 2 4> 
cos 0 


sm </> cos 
s/l - e 2 sin 2 0 


- (1 + e 2 ) [E - (1 - e 2 ) F]J 


1J 

sin 3 

E 

3 1 

„ cos 

E 


s/l - e 2 cos 2 Ef 


(138) 


52 



r b 2 


cos 3 E \/l — e 2 cos 2 E d E = — -jjl + iL(l - e 2 )J [e sin E \/l - e 2 cos J 


E 


+ (1 - e 2 ) In (e sin E +Vl - e 2 cos 2 E)] 


(1 _ e 2 cos 2 E) 3/2 e sin E 
2 



(139) 


/•E, 


s _ , cos 4 E sin E 4 . „ 4 3 _ 

cos 5 E d E = + sm E -- sin 5 E 

5 5 15 


(140) 


f E 2 


E sin E cos 4 E d E = - 


- C ° sS + ll (140) 
^ 5 


(141) 


r 


E sin 3 E d E = — sin 3 E - -I (99) 

2 2 


(142) 


r E 2 


E 2 cos 3 E d E = E 2 I (140) - il (86) 

v 




(143) 


r 


E 2 cos E d E = E 2 sin E - 2 I (86) 


(144) 
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(145) 



E sin 
1 


f E 2 

E 2 cos 3 E d E = I (144) - I (-99) 


|E 2 

E 1 


r E 2 

j E 2 cos E \/l - e 2 cos 2 E d E = I (144) 


— I (145) -^-1 (143) 


(146) 


F 1 

Ecos 5 EdE= — ~ cos 6 E + sin E cos 5, E 

6 36 


+ ® + — sin E cos E + .JL_ sin E cos 3 E 

96 96 144 


(147) 



cos 4 E d E = -E + sin 2 E 
8 16 


cos 3 E sin E 


(148) 


I 


E sin E cos 3 Ed E = - 


E 005 4 E + | I (148) 
4 4 


(149) 
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2 


E 2 cos® E d E = — cos® E sin E + — E 3 
6 48 


+ — E 2 sin £ cos E + E 2 sm E cos 3 E 
16 24 


-I I (92) - ~ I (147) -ll (149) 
o 3 12 



E 2 cos 4 E d E = £2 cos3 E sinE _ & _ i i 


4 4 


( 92 ) 


+ E 2 (E + sin E cos E) - A I (149) 
o 2 



E 2 cos 2 E \/l - e 2 cos 2 E d E = — - I (109) 


-yl (151) - ~I (150) 



cos 5 E d E = 


cos 4 E sin E 
5 



sin E 


4 • 3 r- 

sm J E 

15 




cos 7 £ d E = 


cos 6 E sm E 

7 


+ yl (153) 



(150) 


( 151 ) 


(152) 


( 153 ) 


( 154 ) 
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cos® E d E = COS - 8 - | sin E +(8/9) 1 Cl 54) 
y 


( 155 ) 


E sin E cos 6 E d E = 


■|cos 7 E + il (154) 


>s 6 E dE = I (153) i (154) 


(155) 


E 2 cos 7 E d E = — sin E cos 6 E + — E 2 sin E cos 4 E 
7 35 


+ — E 2 sin E - — E 2 sin 2 E ~ ^ I (141) 
35 35 35 


if I (86) + i| I (142) 


(156) 


E 2 cos 3 E \j 1 - e 2 cos 2 E d E - I (145) 


yl (143-1 -^1 (1585 


e 2 . 



v 4 y 

E sin 5 E d E = --sin 4 EcosE + 1 ^ 

- cos E + ~ cos 3 E 

X 5 5 

L 3 J 


+ J L sinE--sin 3 E + J— sin 5 E 

15 9 25 


e 2 

e i 



r<2 

E 2 sin 4 E cos E d E = — 


si n s E_|I ( 160 ) 


E 


2 


E i 


E 3 sxn 3 E d E = I ( 88 ) - I ( 103 ) 


r E 2 


E 3 sin 3 E cos 2 EdE=^. sm 4 E cos E +-^- 


— cos 3 E - cos E 
3 


-ll ( 161 ) + Li ( 144 ) -II ( 145 ) 
5 5 5 


^ E 2 


E 4 

E 4 sin 2 E cos 3 E d E =— - sm 3 E cos 2 E 

J 


+ 2 - 5 - sin 3 E-ll (15 3 ) - ^ I ( 162 ) 
1 5 5 1 5 



E 4 cos F, d E = E 4 sin E - 4 I ( 88 ) 


(160) 


( 161 ) 


( 162 ) 


( 163 ) 


( 164 ) 


( 165 ) 
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(166) 


E 4 cos 3 E d E = E 4 sin E - 4 I (88) - I (105) 


f 


E 4 cosS E d E = I (166) - I (164) 


E 4 cos E \/l - e 2 cos-* E dE = I (16 5) - — I (166) 


- ~ I (167) 



E 3 sin E cos 3 E d E = - F C ° s4 E + _ I (151) 

4 4 


E 


2 


E 


1 


r E 2 


E 4 cos 4 E d E = j - JL . jE s in E cos E + i E 4 sin E cos 3 E 
40 8 4 


-~I (94) - I (169) 


r e 2 

E 4 cos 2 E d E = gjin E cos E Ef _ 2 j 
1 2 10 


(94) 




E 4 Vl - e 2 cos 2 E d E = 


E s e 2 


— I (171) - T I (170) 


(167) 


(168) 


(169) 


(170) 


(171) 


(172) 
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(173) 


I 


E 3 sin E cos 5 E d E = cos 6 E + \ I (150) 

o 


E 4 cos 6 EdE = — sinE cos 5 E +~E 4 I (148) -~ 
6 o ^ 


.2 i (173 ) 5 1(94 >_!ia69) 

3 4 6 


r^2 


E 4 cos 2 E \/l - cos 2 E d E = I (171) 


^ I (170) ~~ I (174) 
2 ° 


p E 2 

E 3 sin E cos 6 E d E = 


cos 7 E+|l (158) 
7 7 


f 


E 3 sin E cos 4 E d E = I (103) - I (163) 


f** 


E 4 cos 7 E d E = y sin E cos 6 E + A E 4 sin E cos 4 E 


+ 2i^ sinE _^sin 3 E-4l(176) 

35 35 7 


-H® i (88 ) - I (103) I (177) 
35 35 35 


(174) 


(175) 


(176) 


(177) 


(178) 
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62 


E 4 cos 3 E \/l - e 2 co s 2 E d E = I (166) 


I (167) 


I (178) 


cos 8 E d E = 


[cos 7 E sin E 7 
I 7 + T 


cos 5 E sin E 


+ 4-1 (148) 
6 



E 1 

F, sin F cos 8 EdF = — -cos 9 E h — 

9 9 


I (1 55) 


E 


2 


E 


1 


cos 4 E yj 1 - e 2 cos 2 E d E = [l - -^-J I (148) 


5 e 2 \ 


2 4 

— cos 5 E sin E - 5—1 (180) 

12 8 


IE, 


2 


Ecos 7 EsmEdE=- 


E cos 8 E 1 

+ ~I (180) 
8 8 


(179) 


(180) 


(181) 


(182) 


(183) 
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f E 2 

E 2 cos 8 EdE= — sin E cos 7 E + — E 2 sin E cos 5 E 
7 48 


+ E 2 sin E cos 3 E + _^§_E 2 sin E cos E 
192 128 


-14 & --Jl (182) - — I (147) 
384 7 24 


35 

96 


I 


(149) - || I (92) 
64 


(184) 



E 2 cos 4 E \/T 


e 2 cos 3 E d E = I (151) 


-yl (150) - I (183) 



(185) 


/ 2 

-2 9 f j E 2 sin E cos 8 E 8 _ 

I E 2 cos y E d E = + — E 2 s m E cos 6 E 

J 9 63 

E i 


+ E 2 
21 


sin E cos 4 E 4 . _ 4 . 3 _ 

+ - sin E - — sm 3 E 

5 5 IS 


2 

9 


I (181) 


-ra !as61 


128 

105 


I (86) - 


32 

105 


I (141) 


128 

315 


I (142) 


(186) 


61 



E 2 cos 5 E V 1 - e 2 cos 2 E d E = I (143) 


£.1(158) -^.1(186) 

A 8 


(187) 


•r* 


E 3 


E 3 sin E cos 7 E d E = - — cos 8 E 

8 


+ ~I (184) 




(188) 


E 4 cos 8 E d E = sin E cos 7 E E 4 sin E cos 5 E 


48 


+ — sin E cos E + E 4 sin E cos 3 E 
128 128 192 


-^ I( 173 ) -| I( 94 ) -I I(169) 

! 

-il (188) 


f 


E 4 cos 4 E VI - e 2 cos 2 E d E = I (170) 


—■ I (174) - |- I (189) 


(189) 


(190) 


62 



E 3 sin E cos 8 E d E = - 


cos 9 E +ll (186} 
9 3 


**^2 T ?4 q p 4 

E 4 cos 9 E d E = — sin E cos 8 E + —— sin E cos 5 E 
9 63 

T7 


+ M- E 4 sin E cos 4 E + — E 4 sin E 
105 105 


Zll E 4 sm 3 E - — I (17 6} - I (177) 
315 63 105 


25^ I (88) I (162) 

105 315 


-I (191) 
9 


E 4 cos s E VI - e 2 co EdE = I (167) 


— I (178) - ^ I (192) 
2 8 


E 2 cos 2 E d E = — + (— — — ^ sin 2 E 
6 \4 8/ 


E cos 2 E 


(194) 


63 




E 2 VI - e 2 cos 2 EdE =y ~~1 (194) 


4l(151) 




(195) 



sin E (1 - e 2 cos 2 E)" 3/2 d E 


- e cos E 
\J 1 — e 2 cos 2 E 


(196) 



sin E cos E (1 - e 2 cos 2 E)" 3/2 d E = ... 

e Vl - e 2 cos 2 E 


(197) 


sin E cos 2 E (1 - e 2 cos 2 E)~ 3 / 2 d E 




e cos E 


+ sin" 1 ( vl-e 2 cos 2 E) 


%/ 1 ~ e 2 cos 2 E 


■} 


(198) 


IE, 


E 2 sinE (1 - e 2 cos 2 E)~ 3/2 d E = I (87) +-ie 2 I (102) 

% 


+ 2- J I (120) 
8 


E 2 

E 1 


(199) 


E 2 sm E cos E (1 - e 2 cos 2 E) 3/2 d E = I (93) 


,E„ 


+ le 2 I (113) + ^ e 4 I (123) 

^ O 


( 200 ) 


64 



E 2 sm E cos 2 E (1 - e 2 cos 2 Ey 3/2 d E = I (102) 


+ ie 2 I (120) + 11 e 4 I (126) 
2 8 


(201) 


£4 s i n E (1 - e 2 cos 2 E) -3/2 d E = I (107) 

■^1 


+ - e 2 I (104) + — e 4 I (129) 
2 8 


E 2 


( 202 ) 


E 4 sin E cos E (1 - e 2 cos 2 E)" 3/2 d E = I (95) 


+ Ae 2 I (132)+ ile 4 I (131) 

2 8 


(203) 


E 4 sm E cos 2 1 (1 


e 2 cos 2 E) 


f* 3 / 2 d E = I (104) +|-e 2 I (129) 


+ Y e4 1 ( 13 ^ 


(204) 


sin 3 E cos E (1 — e 2 cos 2 E) 3/ ^ 2 d E 


2 - e 2 (1 - cos 2 E) 




\/l - e 2 cos 2 E 


(205) 
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E 2 sin E cos 7 E dE = - 


E 2 cos 8 E l E sin E cos 7 E 
g— + 4 f g 


7 E sin E cos 5 E 35 E sin E cos 3 E 


48 


192 


35 E sin E cos El s ^ 7 „ 

+ z-z-z + — cos* E 4 . ___ cos 0 E 


128 


64 


288 


+ — cos 4 E - I ('91') + — E 2 
768 128 256 


(206) 


f 


E 2 sin 3 E cos 5 E d E = I (123) - I (206) 


(207) 


E 2 sin 3 E cos E (1 - e 2 cos 2 E )~3/2 3 E = I (93) 


- I (113) + Z e2 [I (113) + I (123)] 


+ — e 4 I (207) 
8 


(208) 


66 



r 


V 3 s p J T? E 3 • T 7 7 r> 7 E 3 sin E cos 5 E 

E 3 cos 8 E d E = — sxn E cos 7 E + 

8 48 


+ E 4 + - 35 - E 3 sin E cos 3 E 

512 192 


+ 33 ^ sin E cos E - — I (206 - ) 
128 8 


jgt (123) - ^|l (113) - I (93) [ (209) 


C E 2 


E 4 sin E cos 7 E d E = 


~ — cos 8 E +1 1 (209) 

8 2 


(210) 


E 4 sin 3 E cos E (1 - e 2 cos 2 E) -372 d E = -<I (95) 


I (132) +~e 2 [I (132) - I (131)] 


+ ^e 4 [1 (131) - I (210)] L 


1*2 


8 




( 211 ) 


sin E cos E (1 - e 2 cos 2 E) -1/2 d E = (1 - e 2 cos 2 E) 1/2 


(212) 
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2 


sin 3 E cos E (1 - e 2 cos 2 E) _1/2 d E 


\J 1 - e 2_ cos 2 E 
3~e* 


[(1 - e 2 cos 2 E) - 3 (1 - e 2 )] 


E. 

E 


/* E 2 

E 2 sin E cos 2 E (1 - e 2 cos 2 E)“ 1/2 d E = I (102) 

E 

+ ~I (120) + le 4 I (156) 

2 o 



E 2 sinEcosE(l-e 2 


cos 2 E)" 1/2 dE = I (93) 


+ -Le 2 I (113) *|.e 4 I (123) 
2 8 


r^2 

E 2 sxn 3 E cos E (1 - e 2 cos 2 E)~ 1/2 d E = I (93) 


e 2 

-le 4 I (206) 


(213) 


(214) 


(215) 


(216) 


68 



E 4 sin E cos E (1 - e 2 cos 2 E) d E = I (95) 
Je i 


+ — I (132) + 

2 


3 e 4 


I (131) 


r E 2 


E 4 sin E cos 2 E (1 - e 2 cos 2 E)~^ 2 ‘ d E = I (104) 


+ — I (129) +i-e 4 I (136) 
2 o 


f E 4 sin 3 E cos E (1 - e 2 cos 2 E)" 1 ' 2 d E = I (95) 


+ ^-l)l(132) + (le 4 - e ;)l(131) 
E 2 

-2e 4 I (210) 

8 



sm 3 E cos E (1 - e 2 cos 2 E) 1/f2 


(1 - e 2 cos 2 E) 3//2 


(1 - e 2 cos 2 E) (1 - e 2 ) 


r E 2 


cos 3 E d E = sin E sin 3 E 

3 


(217) 


(218) 


(219) 


( 220 ) 


( 221 ) 
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cos 6 EdE = 


sin E cos 5 E 


£ 

sin E cos 3 E 

24 


+ — sin E cos E 
16 

4 ^ E 
+ T6 

E 2 

E 1 

(*^2 

1“ 

1 

j sin 2 E cos 2 E (1 - e 2 cos 2 E)" 3/2 d E = 

E sin 2 E 
0 + 4 

% 

L 

J 


-|_e* I (153) -ilfll (154) +H?ll (222) 

•*8 8 


- I (221) +|-e 2 I (148) 

4t 


E 2 sxn E (1 - e 2 cos 2 E)' 1/2 d E = I (87) 



2 


I (102) +|-e 4 I (120) 
8 



r 


E 4 sxn E (1 - e 2 cos 2 E)" 1/2 d E = I (107) 


+ ~I (104) +4e 4 I (129) 

2 o 


(222) 


(223) 


(224) 


(225) 


70 



E 2 sin 3 cos E (1 - e 2 cos 2 E) lj/2 dE = I (93) 


-(f* 2 + 1 ) I <U3) -(^-4) I (123, 

4 E2 

+ ~I (206) 


E 4 sin 3 E cos E (1 - e 2 cos 2 E) 1/2 d E = I (95) 


-(y ♦ 1 ) I < 13 «-(T-f) I<13 » 

, e 2 

+ ^ I (210) 

O 


C E 2 


sin E cos 5 E (1 — e 2 cos 2 E)“ 3/2 d E = — 


cos 6 E 


- — e 2 cos® E — — e 4 cos^ E 

16 16 


(226) 


(227) 


(228) 
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I (109) 


E 2 sin 2 E cos 2 E (1 - e 2 cos 2 E)~ 3 / 2 d E — - — 

J 


+ (|e 2 -. ij I (1510 + (~e 4 -le 2 j I (150) 


— e 4 I (184) 
8 


(229) 


E 4 sin 2 E cos 2 E (X - e 2 cos 2 E) -3/2 d E = I (171) 


+ (§e 2 -lj 1(170) + f^e 4 -|e 2 j 1(174)' 


15 


e 4 I (189) 


(230) 


E 2 sin E cos 9 E d E = “-Jq- cos 10 E + 


cos 9 E sin E 


10 


9 /cos 7 E sin E 7 / cos 5 E sin E 5 

10 V 3 8 \ 6 16 


+ JL sm E cos E + — sin E cos 3 E) 

16 24 J 


cos 10 E 9 cos 8 E 7 cos 6 E - 21 E 2 
+ 500 + 3200 + 1600 320 


+ cos 4 E J (91) 

2560 160 


(231) 
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.e 2 


E 2 sin E cos 5 E (1 - e 2 cos 2 E )~ 3/2 dE = ^eM (231) 

8 


+ I (123) +-|-e 2 I (206) 


(232) 


E 4 sin E cos 9 E d F = - 11 cos 10 E + | * _ 0 -° .. i l s:Ln E 


10 


cos 7 E sin E 7 /cos 5 E sin E 5 ... 

8 8 l 6 16 


5 5 \ 

sin E cos E + — s in E cos 3 E 

16 24 ' 


-^i(206) -^icu3) 


189 £2 - i 89 -! (93) - — I (113) 
1280 640 320 


(233) 


/• E 2 


E 4 sin E cos 5 E (1 - e 2 cos 2 E)" 3/2 d E = — e 4 I (233) 

8 


+ I (131) -E-ye 2 I (210) 


(234) 
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sin 4 (1 - e 2 cos 2 E) 1/2 , ^ 1 . , _ 

1 dE = - — sin 3 E- sin E 

cos E 3 


- i_ s i n S E _ 

10 8 


i sin 3 E (— + — cos : 
7 l5 5 


- cos 4 E 


'E, 


f 2 ■ 2 

si,n - £ Cl - e 2 cos 2 E) 1/2 d E = - sin E-ie 2 I (97) 
„ r cos E 2 


ism 3 E (1 + cos 2 E 
5 V 3 


^2 


E sin 7 EdEi. 


~e 


E sin 4 E cos E / . „ 6 ' 

/sin 2 E + — 


+ ~ E cos E (cos 2 E - 3) +i ( 1 sin 7 E 


+ ~ sin 5 E - JL sm 3 E + ^ sin Ej 


r 


E 2 sin* e cos 3 EdE=il (237) - ~ sin 7 E 


+ I (161) 


E 


(235) 


(236) 


(237) 


(238) 
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2 


£2 s i - n4 1 (1 - e 2 cos 2 E) 1 ' 2 d E = - -f! I (238) 
cos E 8 

e 2 

( 239 ) 

E i 


-1 e 2 I (161) - I (144) - I (99) 


f 2 

E 4 sin 4 E cos 3 E d E = I (164) - I (167) 

% 


+ I (178) 


(240) 


r 


E 4 sin 4 E cos E d E = I (105) - I (164) 


(241) 


r 


E 4 sin 4 E 
cos E 


d E= I (165) - I (105) 


(242) 




E 2 sin 4 E (1 - e 2 cos 2 E) 1/2 
cos E 


d E = I (242) 


- I (241) -£ll (240) 
Z o 


(243) 
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E 2 sin 2 E (1 - e 2 cos 2 E) 1/2 
cos E 


d E = - — I (99) 
2 


- I (144) -i-I (119) 
o 




E 4 sm 2 E (1 - e 2 cos 2 E) 1 / 2 
cos E 


I (105) 


| I (164) - I (165) 


(244) 


(245) 


Thus, both terms of equation (6. a) are known for each of the Vinti orbital ele- 
ments, and the air-drag effect for any eccentricity is also known. 

REMARKS 

The derivation of the drag rate equations is perfectly general and stems _ 
only from first principles. That is, since we employ an angular momentum 
argument, these equations represent a time deviation of the orbital elements 
from forces conservative or nonconservative. Here we have solved them from 
semiempirically determined force functions of atmospheric drag through the 
eccentric anomaly, hi all of our test cases, we used only a static model 
atmosphere. Next is the inclusion of a more refined day- night profile. 

The integrated rate equations now permit one to include drag by reinitiali- 
zation; after solving Yinti's kinetic equations of motion for a given value of time, 
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the computed orbital parameters are then used to obtain the solutions of the rate 
equations for that same time. The latter results are then used to obtain the new 
constants of the motion in Vinti’s theory which arise in the factorization of cer- 
tain quartic polynomials, necessary in the separation of the Hamilton- Jacobi 
equation. This entire process is carried out for each given time. It is con- 
firmed that for satellites encountering sufficient aerodynamic drag, this pro- 
cedure adds considerable accuracy in the drag region, and allows one to perform 
fewer corrections on the orbit. It must be emphasized, that this procedure is 
entirely different from that of numerical integration. In the Vinti program with 
aerodynamic drag, the rate equations serve only to provide new initial condi- 
tions as functions of time through the orbital elements, with which separability 
can be achieved, and the kinetic equations of motion can be solved (Reference 2). 

A Fortran version of this program has been prepared and tested on the 
IBM 360/95 electronic digital computer. At present, the Fortran program does 
not include the second term of equation (6. a), (Equations (86) through (245)) which 
account for drag variations on the satellite from changing atmospheric density. 
This part is presently being added to the Vinti orbit determination program. On 
the other hand, if the eccentricity is not too high, the second term of equation 
(6, a) may not have a great effect, providing perigee is not low also. In any 
event, it is interesting to observe the effect of even the first term in equation 
(6. a) on the accuracy of a heavy air-drag satellite. Of concern are the San 
Marco 2 and S-55 Micrometeor Impact Satellites. To date, only the San Marco 2 
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satellite has been tested. The initial orbital data for this satellite are as 


follows: 


Semimajor axis 
Eccentricity 
Perigee 
Apogee 


6807.0374 km 
0.03333142 
201.7898 km 
655.5663 km 


Inclination 2.89095642 degrees 

Period 93.15 minutes 

Figure 1 is a plot of the time variation of the Izsak elements after reinitializing 
the Vinti program using these initial conditions. Their behavior is exactly what 
„ is expected (see Sterne, Reference 1). Figures 2 and 3 are the size of residuals 
of computed and observed direction cosine differences AL and A M versus time. 
The comparison was for a three day period. In both cases the non- drag and drag 
orbits were comparable for approximately twenty hours. After this, due to heavy 
air drag, the difference was large, as much as an order of magnitude. En the drag 
calculations, only a static atmospheric model was used. However, this is being 
replaced by a more refined day-night profile model with which somewhat more 
accuracy should be obtained. At any rate, the results are striking, and it is felt 
that with the addition of equations (86) through (245) (the second term of equation 
(6 .a)), there exists an analytic drag-gravitational theory or program operating 
as a fundamental umt to account for the motion of an artificial Earth satellite. 

The drag program has also been incorporated into the differential correc- 
tion portion of the Vinti program. This has added the following advantages: Since 
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from Figures 2 and 3 it is seen that the Vinti orbit generator can predict 
approximately three and one-half times longer with the drag added, an orbit 
correction to obtain mean elements can be made using the seventy-two hour arc 
instead of the twenty hour arc. Secondly, the mean fitted elements now reflect 
or contain the characteristics of atmospheric perturbations as well as those of 
the full Vinti gravitational potential. There results a more accurate prediction 
over a much longer period of time, without the need of frequent orbit correc- 
tions. This program is presently being prepared for testing the S-55 Micro- 
meteor Impact "satellite. 
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APPENDIX 


As mentioned in Section n above, we shall give here some examples on the 
solution of the types of integrals encountered there. Basically, all of the differ' 
ent integrals are approached in the same manner as we shall show here. 

Given the integral, 


Let 



cos 2 E) 1/2 dE 


x 2 - 1 - e 2 cos 2 E 


dE 


x dx 

e 2 cos E sinE 


cos 2 E 



sin 2 E 


e 2 - I +x 2 


e 


2 


let 


k 2 = 
sin 2 E 


1 - e 2 
x 2 — k 2 


81 




(1 


~„2 


cos 2 E) 1/2 dE 



x 2 dx 

(X -x 2 ) 1/2 (x 2 -k 2 ) 


E(a, 4) 


from Reference 3. For the integral, 



sm 3 E cos E( 1 + e cos E) ^ 2 


( 1 - e cos E)~ 3/2 dE 


■ r 


sin 3 E cos E(1 + e cos E) 

( 1 - e cos E) (1 ~ e cos E) 


1/2 


1/2 


dE 


multiply numerator and denominator by ( 1 - e cos E) 1/2 



sin 3 E cos E(1 + e cos E) 1/2 (1 - e cos E) 1/2 
( 1 - e cos E) ( 1 - e cos E)^ 2 (1 - e cos E) 1/2 


■f 


2 sm 3 E cos E( 1 - e 2 cos 2 E) 
(1 - e cos E) 2 


1/2 


dE 
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Now use the same substitutions of example 1 (x 2 = 1 - e 2 cos 2 E) 


f 


sin 3 E-cos Ex 2 dx 


x e 2 sxn E cos E [l ~ ( 1 - x 2 ) 1/2 ] 


1 f 2 sin 2 j 

e2 L fi-ci- 


sm 2 E x 2 dx 




x 2 -k 2 

•E. — * X 2 dK 


j_ r 2 e 

e2 L fl- 


! [l “(1 -X 2 ) 1/2 ]‘ 


1 f 2 (x 2 - 1 

e4 Je, fl-a- 


- k 2 ) x 2 


, [l-(l-x 2 )'«]' 


dx 



f 2 x 4 dx 

k 2 1 

p 2 X 2 dx 

e 4 . 

E x [l -(1 ~X 2 ) 1/2 ] 

| 2 ~ e 4 j 

L, 


multiply numerator and denominator by 


[l + ( 1 - x 2 



2 
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* E 2 

x 4 | 

1 + (1 

- X 

.2)1/2’ 

2 

dx 

e 4 . 

e i 1 

[l-(l -> 

r 2) 1/2] 

2 | 

[l + (1 

_ x 2) 1 /2] 


kf_ 

* E 2 

X 2 

1 +(1 

_ x 2) 1/ 2’ 

2 

dx 

e 4 

E 1 1 


■ 2 ) 1/2 ] 

| 2 | 

[1 +(1 

-X 2 ) 


1/2 


Now use the fact that (a - b) 2 (a + b) 2 = (a 2 -b 2 ) 2 


_1_ 

r E 2 X 4 

1 

+ (l-x 2 ) 1/2 

1 1 

e 4 . 

E 1 1 

[l-(l -X 2 )] 

2 


Jk^ 

" E 2 x 2 [ 

1 + (l-x 2 ) 1/2 

] 2 dx 

e 4 . 

E i 

[l -(1 -X 2 )] 

2 


1 

' E 2 X 4 

1 + 2( 1 - x 2 ) 1/2 + ( 1 -x 2 ) j 

dx 

e 4 . 

E 1 

x 4 



k 2 

' E 2 X 2 

1+ 2( 1 - x 2 ) 1/2 + (1 - x 2 )] 

dx 

e 4 . 

E i 

X 4 




[2 +2(1 -x 2 ) 1/2 - x 2 ] dx 


k 2 r E 2 1 

~~ ~ [2 +2(1 -x 2 ) 172 ~ x 2 ] dx 

e x" 1 L 
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2(1 ~x 2 ) 1/2 dx 


,e 2 

| x 2 dx - k 2 
Je i 




>*E f 


•'Ei 


r*' 


x 2 dx - 2k 2 


dx 

T 2 


- 2k 2 



(1 — x 2 ) 1/2 dx 
x 2 ( 


“7 H 2 +k2 ) x + 2 • T [ x (! - 


X 


2 s 1/2 


+ sin 




X~ 

~T + 2k^ 


+ 2k 2 


(1 -x 2 ) 172 
+ sin x 

_ 
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Now substitute x 2 - 1 -e 2 cos 2 £ back in for x 
= ~ |(1 -e 2 cos 2 E) i/2 (3 -e 2 + e cos E) 

+ (3 - 2e 2 ) sin" 1 [(1 -e 2 cos 2 E) 1/2 ] 

- 3 -(l-e 2 c°s 2 E) 3/2 + 2 (l- e2) 

1 E 2 

x (1 -e 2 cos 2 E)" 1/2 (1 +e cosE)> ^ 
sm 3 E cos E( 1 + e cos E) i/2 ( 1 - e cos E)" 3/2 dE . 

Finally, for 

f E 2 

e cos E( 1 - e 2 cos 2 E) 1/2 dE 
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use the same substitution as used in first two examples. Then this becomes 


C*” 2 e cos Ex 2 dx 

f E2 X 2 dx 

f 2 x 2 dx 

J e 2 cos E sinE 

E i 

J P e(x 2 -k 2 ) 1/2 

1 e 

J Ej <K=-k2)‘« 


Now let 


x = k sec t 


dx - k sec t tan t dt 



k 2 sec 2 t k sec t tan t dt 
[k 2 sec 2 t - k 2 ] 1/2 


k 3 sec 3 t tan t dt 
k [sec 2 t - l] 



k 3 sec 3 t tan t dt 
k( tan 2 t ) 1/2 


k 2 sec 3 t dt 


k 2 [sec t tan t + ln(sec t + tan t)]^- 
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Y /2 i 

- 1 = Y (* 2 -k 2 ) 1/2 



Substitute back in for 


x 2 - 1 - e 2 cos 2 E 


x 

1 tan t 


(sec 2 t - 1') 1/2 = f — 


88 



2 


|e smE(l 


- e 2 


cos 2 E) 


1/2 


+ 


( 1 ~ e 2 ) In 


sinE+(l~e 2 cos 2 E) 
' ( 1 - e 2 ) 1 / 2 




e cos E(1 - e 2 


cos 2 E) 1//2 dE . 


These three examples give the general methods used in solving the 
integrals. 
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Figure 1. San Marco 2 Satellite. Variation of Elements 
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Figure 2. San Marco 2 Satellite (AL vs. T) 










